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Research Objectives 

A quantitative understanding of photolysis rate coefficients (or “y-values”) is essential to deter- 
mining the photochemical reaction rates that define ozone loss and other crucial processes in the 
atmosphere, j - Values can be calculated with radiative transfer models, derived from actinic flux 
observations, or inferred from trace gas measurements. The principal objective of this study is to 
cross-validate /-values from the Composition and Photodissociative Flux Measurement (CPFM) 
instrument during the Photochemistry of Ozone Loss in the Arctic Region In Summer (POLARIS) 
and SAGE III Ozone Loss and Validation Experiment (SOLVE) field campaigns with model cal- 
culations and other measurements and to use this detailed analysis to improve our ability to de- 
termine /-values. Another objective is to analyze the spectral flux from the CPFM (not just the 
j- values) and, using a multi-wavelength/multi-species spectral fitting technique, determine atmos- 
pheric composition. 


Progress and Results 

1. Improved Agreement of Modeled and Measured j- values in the Lower Stratosphere 

This study begins with detailed comparisons of modeling and measurements made during the In- 
ternational Photolysis Frequency Measurement and Modeling Intercomparison (IPMMI). IPMMI 
was conducted during a one-week intensive at the Marshall field site of the National Center for 
Atmospheric Research (NCAR), near Boulder, Colorado, June 15-19, 1998, followed by analysis 
over several years [ Cantrell et al., 2003]. Twenty-one researchers from eight different institutions 
participated in the measurement portion of the intercomparison, and 18 modelers from 13 groups 
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Table 1 . N0 2 cross sections and extraterrestrial solar flux used in initial and revised versions of 
the data 


y- Value Source 

N0 2 Cross Section 

ET Flux 

APL radiative transfer model 

initial (POLARIS version) 

DeMore et al. [1997] 

MODTRAN 

revised 

Harder et al. [1997] 

ATLAS 

IPMMI spectroradiometer (/no,) 

initial ( DeMore et al.) 

DeMore et al. [1997] 

— 

revised ( Harder et al.) 

Harder et al. [1997] 

— 

CPFM spectroradiometer (/no,) 

initial (POLARIS) 

DeMore et al. [1997] 

— 

revised (estimated) 

Harder et al. [1997] 

— 


provided modeling calculations. The intercomparison focused on spectral actinic flux [ Bais et al., 
2003] and the photolysis rate coefficients for N0 2 , 0 [Shetter et al., 2003], and ozone photo- 

lysis to O('D), 7o,->0(‘D) [Hofzumahaus et al., 2004]. One of the major goals of IPMMI was 
to make recommendations as to which of the various model input data (e.g., extraterrestrial solar 
irradiance, absorption cross sections) lead to the greatest consistency between measurements and 
modeling calculations, based on the results of the intercomparison. Analysis was supported in 
part by this grant because of the significant impact the IPMMI findings have on the calculation 
of /-values using CPFM data and the comparison of CPFM /-values with both models and in situ 
measurements of NO* partitioning. 

The results described in the following paragraphs are summarized in Tables 1 and 2. 

Extraterrestrial solar flux. A number of models in IPMMI used the ATLAS extraterrestrial 
solar spectrum [Kaye and Miller, 1996; Woods et al., 1996]. The APL model, in contrast, used the 
MODTRAN3 spectrum [Berk et al., 1989], as was used for the APL model-CPFM comparisons 
reported by Swartz et al. [1999] and Del Negro et al. [1999] during the POLARIS campaign. 
This difference led to some significant deviations from the other models and with actinic flux 
measurements at the ground. In particular, over 300-315 nm, where ozone photolysis occurs, 
the MODTRAN flux exceeds ATLAS by roughly 7%, with excursions of up to 20%. Thus, the 
selection of the solar spectrum alone increased APL modeled flux by several percent in this region, 
which is important for ozone photolysis. 

The ATLAS-derived model calculations agreed better with the measured spectra at the ground 
during IPMMI, which suggests the superiority of the ATLAS spectrum for model calculations 
[Swartz, 2002; Bais et al, 2003], This also supports the recent general trend found in the literature 
toward using the ATLAS extraterrestrial flux [Grobner and Kerr, 2001]. The ATLAS flux (290- 
420 nm) should be used in y'-value modeling. 

Support for larger N0 2 cross sections. All of the models used in IPMMI (including the APL 
model) utilized the DeMore et al. [1997] or comparable N0 2 cross sections, (Jncv The models 
were compared to j NOi chemical actinometers in order to determine the absolute accuracy of the 
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size of the N0 2 cross sections. If a rough uncertainty of 10% for the model calculations is as- 
sumed [Shetter et al, 1992; Weihs and Webb, 1997], along with 7% for the actinometers, only nine 
of the 15 models used in IPMMI agree with the measured 7 NC , 2 to within the combined model- 
measurement uncertainty [Swartz, 2002], If larger N0 2 cross sections were used, however, like 
those of Harder et al. [1997], which lead to a j N o 2 of about 7% higher, all of the models would 
agree with the measurements to within the combined uncertainties, and eight of the models would 
agree to within the measurement uncertainty alone. This finding is consistent with the results of 
better IPMMI spectroradiometer-actinometer agreement when using the Harder et al. [1997] cross 
sections [ Swartz , 2002; Shetter et al., 2003]. 

This conclusion also confirms results from the limited spectroradiometer-actinometer compar- 
isons of Kraus et al. [2000], which showed better agreement when using the larger Merienne et al. 
[1995] cross sections, in comparison with those of DeMore et al. [1997]. The community should 
thus consider using larger N0 2 cross sections [e.g., Merienne et al., 1995; Harder et al., 1997] 
and/or further laboratory evaluations of o N q 2 . 

Implications of IPMMI findings for CPFM and POLARIS data. The results from the IPMMI 
intercomparison have important implications for the CPFM. Replacing the MODTRAN extrater- 
restrial flux with the ATLAS data reduces APL model j 0 ^ by about 7% on average. This impressive 
correction brought the APL model into nearly perfect average agreement with y 0 , chemical acti- 
nometer and spectroradiometer measurements during IPMMI. If the POLARIS j Q ^ calculations 
[Swartz et al., 1999] were re-computed with the ATLAS extraterrestrial flux, the APL CPFM 1 /CPFM 
agreement would go from +14% to about +7%, and the APL TOMS 1 /CPFM agreement would drop 
from +8% to +1%, dramatically improving agreement with CPFM flux-derived j- values [Swartz, 
2002 ], 

Since both the APL model and CPFM used the same cross section for j NOn during POLARIS 
[Swartz et al., 1999], increasing cr NOi would not affect their agreement. However, Del Negro et al. 
[1999] compared modeled and CPFM-derived j n02 to steady-state y N o 2 , 7 no 2 ,ss’ inferred from 
measured NO, N0 2 , 0 3 , CIO, and H0 2 . Use of the larger Harder et al. [1997] N0 2 cross sections 
during POLARIS would increase the model and CPFM-derived y N o 2 , improving agreement with 
7no 2 ,ss- The same is true for [N0 2 ] inferred from in situ trace gas measurements and modeled 
and CPFM — using larger N0 2 cross sections would significantly improve CPFM and APL 
agreement with the measured chemical composition along the ER-2 flight track. 

We are currently preparing a manuscript for publication (Swartz, W. H., R. E. Shetter, L. A. 
Del Negro, and C. T. McElroy, “Improved agreement of modeled and measured y'-values in the 
lower stratosphere,” in preparation), which will include our reevaluation of the POLARIS j- values, 
modeled, derived from CPFM, and inferred from in situ chemistry Ono 2 ,ss from Del Negro et al. 
[1999]). 

The high-solar zenith angle conditions during SOLVE posed a much greater challenge. Un- 
fortunately, technical difficulties rendered the CPFM data from SOLVE ultimately not sufficiently 
tractable to warrant further detailed analysis of the CPFM j-value algorithm in the context of 
SOLVE. We thus turned our attention instead to the analysis of CPFM spectral flux data from 
POLARIS to support the development of a multi-wavelength/multi-species spectral fitting retrieval, 

1 APL model calculations based on CPFM and TOMS column ozone and surface reflectivity data, APL cpfm and 
APL toms , are described by Swartz et al. [1999]. 
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which would be of tremendous importance in analyzing datasets from new aircraft- and space- 
based platforms (e.g., OMI). 

2. Multi-Spectral Retrieval of Ozone Column using CPFM Data 

Overview. The CPFM spectrometer measures the direct attenuated solar irradiance across a hor- 
izontal diffuser plate on the top of the instrument [McElroy, 1 995], The density of the atmosphere 
and resultant scattering above ER-2 altitudes is very small. Thus, the measured irradiance, when 
cosine-corrected for the solar zenith angle, is essentially equal to the direct solar irradiance — 
assumed the total downwelling portion the radiation field. In this study we have used the upward- 
looking CPFM data to infer the ozone column between the instrument and the Sun. We have 
used standard least squares fitting techniques to perform this inversion. In order to do so, we con- 
structed a forward model of the observations — that is, a functional relationship between the CPFM 
measurements and various free parameters, including the ozone column. Standard fitting routines 
were then used to determine the values of the free parameters most likely to produce the CPFM 
measurements. 

Our ability to infer the ozone column is limited by the fidelity of the forward model. The 
model fidelity, in turn, is limited by our knowledge of the instrument, as well as various physical 
quantities such as the solar spectrum and the extinction cross sections involved. Our knowledge 
of the instrument is based on our collaboration with C. T. McElroy (Environment Canada). Our 
model of solar occultation is based on the stellar occultation work described by Yee et al. [2002] 
and DeMajistre and Yee [2002], 

We will first describe the forward model, including both the occultation and instrument char- 
acterization. The retrieval method will then be explained. 

Forward model. Each spectral element of the CPFM data, S -, can be written as 

Sj = Aj( u) J s(X - \)I(X)dk, (1) 

where /(A) is the directly transmitted solar radiance entering the CPFM instrument, .v(A — A ( ) is the 
instrument spectral spread function, and A -(u) is the CPFM angular response, i.e., the sensitivity of 
the instrument for a given unit vector to the Sun, u, measured in instrument coordinates. The solar 
radiance is independent of the instrument and is treated separately from the instrument functions 
Aj( u) and s(A — A,) below. 

Radiance model . The solar radiance entering the CPFM can be written as 

/»oo 

/(A) = /„(A)exp(-^ / <y k n k ds), (2) 

k Jo 

where /^(A) is the extraterrestrial solar irradiance, and a k and n k are the extinction cross sections 
and corresponding densities for the constituents along the path to the Sun, ds. The constituents 
considered are (1) 0 2 and N 2 (providing Rayleigh scattering), (2) ozone, (3) aerosols, (4) N0 2 , 
and (5) N0 3 . At some wavelengths for some constituents (including ozone), the cross section 
depends on the temperature and is therefore, strictly speaking, a function of position. We assume 
that the temperature measured on the aircraft is representative of the temperature for most of the 
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absorbing constituents along the path (this assumption will result in a small systematic error). We 
now write 

I(X) = L.Wa f (-J,aiN t ), (3) 

k 

where o£ are the effective cross sections and N k are the column densities for the constituents k. 
This is the forward model that we use for the radiance at the CPFM. 

CPFM model . The instrument spectral spread function, s ( A - A ( ), was provided by C. T. McEl- 
roy and has been used directly. The treatment of the spatial response, A -(u), is a bit more compli- 
cated. The intended design of the CPFM was to have the spatial response be proportional to cos 0, 
where 0 is the solar zenith angle when the instrument is in its nominal orientation. To the extent to 
which this goal is met, the CPFM measurements are proportional to the downwelling hemispheric 
flux. For this reason, the CPFM data have been reported as 



Cj cos 0 ’ 


(4) 


where F } is the downwelling flux and Cj is a calibration constant. To a small extent, however, 
this design goal was not met, particularly at longer wavelengths. The discrepancy is a rather 
complicated function of CPFM orientation for which a general characterization is not available. 
However, Tom McElroy has provided us with an estimate of the radiance for the September 2 1 , 
1997, flight where this effect has been accounted for. We have used this data set to calculate a 
revised hemispheric flux, F* , 


= fs{k-\)I{X)dX, (5) 

J cos 0 7 

which are the data used in the retrievals. 


Retrieval. Combining (3) and (5), we can write 

F; = -L- [ s(X - \)U{X)e-^k N kdX. (6) 

J COS v J 

We use a modified version of this equation in the retrieval in which we define the free parameters 
explicitly: 


F* = [s(X - UUA)e- {a ' N °K* y^y+^o^+^o^o+^^dX, (7) 

J cos 0 7 ' 

where the N° are climatological values of the column densities and the a n are the free parame- 
ters to be retrieved. The summation in this equation employs climatological values of the other 
constituents along the path. The parameter a {) when used in the retrieval accounts for forward 
model errors in calibration, aerosol extinction, and other factors that have a proportional effect on 
the corrected hemispheric flux. This final equation (7) is used in a standard Levenberg-Marquardt 
non-linear least squares fitting procedure. The results are reported below. 
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Results. The ER-2 transit flight from Fairbanks, Alaska to Barber’s Point, Hawaii on September 
21, 1997, was used for this analysis. Figure la shows the solar zenith angle and aircraft altitude 
along the flight track. In Figure lb, the CPFM spectrum measured at the minimum solar zenith 
angle for the flight (35.6°), corrected for the “non-ideal” response of the diffuser plate, is shown, 
along with the best fit from the APL model. From this best fit the column ozone was inferred. 

Figure 2 shows the retrieved ozone column throughout the flight. For comparison, the values 
retrieved by McElroy (C. T. McElroy, unpublished data) are also shown, along with the predicted 
ozone column. The ozone column was predicted assuming a horizontally homogeneous ozone 
field (profile based on ozonesonde data and scaled to the TOMS total ozone along the flight track 
[Swartz et al., 1999]) and accounting for atmospheric refraction (based on 3-dimensional NCEP 
reanalysis data [ Kalnay et al., 1996]). The overall agreement is excellent, with the APL retrieval 
agreeing with our predicted column to within roughly 5%, except for the beginning of the flight. A 
number of factors may be causing discrepancies, including uncharacterized angular response of the 
CPFM diffuser and inaccuracies in the assumed ozone field (used for the predicted ozone column, 
including spatial inhomogeneities). 

This multi- wavelength/ multi-species spectral fitting retrieval based on CPFM data represents a 
preliminary effort. We are currently seeking opportunities to continue this work, analyzing CPFM 
and other multi-spectral datasets. 


Summary 

/- Value comparisons among models and measurements during IPMMI and POLARIS demonstrate 
the superiority of the ATLAS extraterrestrial solar spectrum (290-420 nm), which should be used 
for /-value modeling. Comparisons from IPMMI and POLARIS also show better agreement when 
using larger N0 2 cross sections [e.g., Harder et al., 1997]. The community should thus consider 
using larger N0 2 cross sections and/or further laboratory evaluation of ct n0 . The use of the 
ATLAS extraterrestrial flux and larger N0 2 cross sections would remove a significant portion of 
the discrepancies of CPFM /-values with both the APL model and with measured NO A partitioning 
during POLARIS (see initial and revised comparisons in Table 2). This also suggests that the 
algorithm currently used to convolve CPFM flux measurements and derive /-values is adequate 
for /-value calculation under POLARIS ER-2 conditions. We are currently preparing a manuscript 
for publication (Swartz, W. H., R. E. Shelter, L. A. Del Negro, and C. T. McElroy, “Improved 
agreement of modeled and measured /-values in the lower stratosphere,” in preparation), which 
will include our reevaluation of the POLARIS /-values modeled, derived from CPFM, and inferred 
from in situ chemistry (Jno 2 ,ss f rom D e l Negro et al. [1999]). 

A preliminary multi-wavelength/multi-species spectral fitting retrieval analysis of the CPFM 
POLARIS data shows great promise, and we hope to apply it to other datasets (such as OMI) in 
the near future. 
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Time (sec UT) 



Wavelength (nm) 

Figure 1 . (a) Flight solar zenith angle (SZA) andER-2 altitude for the September 21, 1997, transit 
flight from Fairbanks, Alaska to Barber’s Point, Hawaii. The vertical dotted line represents the 
point of smallest SZA (35.6°), at 84,040 sec UT. (b) Measured and model fitted solar irradiance 
spectra at 84,040 sec UT (SZA=35.6°). The CPFM measurements have been corrected for the 
“non-ideal” behavior of the CPFM diffuser plate, but not scaled by sec(SZA). The best APL 
model fit of the spectrum is overlaid. 
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Time (sec UT) 



7.0x1 0 4 7.5x1 0 4 8.0x1 0 4 8.5x10 4 9.0x10 4 

Time (sec UT) 

Figure 2. (a) Predicted and retrieved slant ozone column using the APL and McElroy re- 
trievals. (b) The fractional difference between the two retrieval techniques and calculations: 
(retrieval/predicted) — 1. 
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